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Abstract In these lectures we describe how a theory of quantum gravity may be 
^f) ' constructed in terms of a lattice formulation based on so-called causal dynamical 

triangulations (CDT). We discuss how the continuum limit can be obtained and 
how to define and measure diffeomorphism-invariant correlators. In four dimen- 
sions, which has our main interest, the lattice theory has an infrared limit which 
can be identified with de Sitter spacetime. We explain why this infrared property 
of the quantum spacetime is nontrivial and due to "entropic" effects encoded in the 
pH ' nonperturbative path integral measure. This makes the appearance of the de Sitter 

universe an example of true emergence of classicality from microscopic quantum 

laws. We also discuss nontrivial aspects of the UV behaviour, and show how to in- 

^ , vestigate quantum fluctuations around the emergent background geometry. Finally, 

f^i ' we consider the connection to the asymptotic safety scenario, and derive from it a 

VO , new, conjectured scaling relation in CDT quantum gravity. 
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1 Introduction 

How to reconcile the classical theory of general relativity with quantum mechanics 
remains an unsolved problem. Flat Minkowskian spacetime seems an excellent local 
approximation to spacetime down to the smallest distances we can probe in the lab- 
oratory. At least the Standard Model of elementary particles, which relies heavily on 
the Minkowskian spacetime structure, works almost too well, leaving us presently 
with little clue as to what should replace it at shorter spacetime distances, i.e. higher 
energies. If one naively tries to quantize the theory of general relativity by making a 
perturbative expansion around this flat background one is faced with the fact that the 
corresponding field theory is not renormalizable. The mass dimension of the gravi- 
tational coupling constant is —2 in units where c = h = 1. To deal with this problem 
one can try to go beyond conventional quantum field theory. One such attempt is 
string theory. However, until now it has added little to our understanding of why 
to a very good approximation we live in a 3+1 dimensional classical world gov- 
erned by Einstein's equations with a positive cosmological constant, around which 
there presumably are small quantum fluctuations. Loop quantum gravity is another 
attempt to quantize gravity, which introduces new ways of treating gravity at the 
Planck scale, but it has problems with recovering classical gravity in the infrared 
limit. Here we will describe a much more mundane approach using only standard 
quantum field theory. In a sum-over-histories approach we will attempt to define a 
nonperturbative quantum field theory which has as its infrared limit ordinary clas- 
sical general relativity and at the same time has a nontrivial ultraviolet limit. From 
this point of view it is close in spirit to the renormalization group approach, whose 
application to gravity with the hope of establishing its asymptotic safety was first 
advocated long ago by Weinberg [1], and more recently substantiated by several 
groups of researchers [2]. 

The approach reported here is nontrivial for two reasons which combine to make 
it genuinely nonperturbative. First of all, as just stated, a theory of quantum gravity 
is not perturbatively renormalizable, and thus, whatever field theory one invents, it 
must in some sense be nonperturbative. Here we want to use a lattice to provide an 
ultraviolet regularization of the quantum field theory. The lattice regularization of 
quantum field theories has been very successful, but it is usually implemented in flat, 
Euclidean spacetime, where the Osterwalder-Schrader axioms ensure us that this is 
unproblematic. One knows how to get from a quantum field theory formulated in 
Euclidean spacetime to a quantum field theory formulated in Minkowskian space- 
time. However, little is known about the analogous issue once we move from flat to 
curved spacetimes and even more, to a situation where spacetime itself becomes the 
object of quantization 1 . That the situation is nontrivial even at the most elementary 
level is seen by considering the Einstein-Hilbert action from which we can derive 
the classical equations of general relativity. It is formally straightforward to rotate 
this action from Minkowskian to Euclidean signature. However, one then has to face 



1 An extension of the Osterwald-Schrader axioms to certain diffeomorphism-invariant theories was 
given in [3]. 
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the fact that the Euclidean action is unbounded below, the unboundedness caused by 
the "wrong sign" of the conformal mode, corresponding to overall, local rescalings 
of the metric tensor. This is a potential obstacle for constructing the quantum theory: 
when trying to sum over all geometries, with weights given by the exponential of 
minus the action, there are geometries with arbitrarily large negative action which 
may render the sum over paths ill defined and the Euclidean theory nonsensical. 
Consequently, the UV lattice regularization of the path integral has to be such that 
it also regularizes the infinities which can arise due to the conformal factor. Even if 
such a regularization exists (and it does, as we shall see), how can one ever expect to 
obtain something finite in the continuum limit where the regularization is supposed 
to be removed? 

The only solution from a continuum point of view is that the correct path integral 
measure in the Euclidean sector suppresses the unbounded conformal factor. In the 
lattice approach the measure factor is of an "entropic" nature: it reflects how many 
configurations (microscopic, geometrical realizations) there are corresponding to a 
given value of the action. This entropic factor will enter as an integral part of the 
bare effective lattice action. (We will illustrate this below in toy examples where 
everything can be calculated analytically.) The "entropy part" of the effective action 
will be independent of its "bare coupling constant part". Usually, the possibility of 
obtaining a continuum limit of a lattice theory is linked to the existence of a critical 
point (more generally, a critical surface) in coupling constant space. This will also 
be the case here. The entropy part of the effective action plays a crucial role in deter- 
mining the critical value of the bare coupling constants, and the continuum quantum 
field theory will then emerge at that critical point. However, in such cases there may 
be no "obvious" continuum theory one can read off from the lattice effective action 
since we might not know the precise form of the entropy part of the effective action. 
This highlights the truly nonperturbative nature of the continuum theory. 

This review article' is organized as follows: first we describe how a two- 
dimensional toy model of quantum gravity can be solved explicitly, illustrating some 
of the points made above. Then we describe the four-dimensional theory, the nu- 
merical results obtained by Monte Carlo simulations and how to connect the lattice 
formalism to the renormalization group approach and to a new theory, so-called 
Hof ava-Lifshitz gravity. 



2 The CDT formalism in two dimensions 
2.1 Generalities 

The lattice formulation of Euclidean quantum gravity, i.e. the quantum theory of 
Euclidean geometries, has been very successful in two dimensions. In two dimen- 
sions, gravity does not have any field-theoretic degrees of freedom, but neverthe- 



2 For less technical accounts of this approach to quantum gravity, see [37, 38]. 
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less two-dimensional Euclidean quantum gravity constitutes a nontrivial example 
of a diffeomorphism-invariant quantum theory of geometries. The lattice theory, 
regularized by the method of so-called dynamical triangulations (DT), provides a 
diffeomorphism-invariant cut-off of two-dimensional Euclidean quantum gravity. It 
is thus a misconception that a lattice regularization will necessarily break diffeo- 
morphism invariance. Rather, one should view the use of DT in the path integral as 
a way to sum directly over geometries, thus avoiding completely the issue of dif- 
feomorphism invariance. The reason why such an interpretation is possible is that 
the triangulations used in DT can be viewed as piecewise linear geometries without 
any specific metric assigned to them: once we know the lengths of the links and the 
gluing of the simplices, we have the complete information about the geometry. Us- 
ing identical simplices the basic information about the geometry is entirely encoded 
in the way the simplices are glued together, and the summation over geometries be- 
comes the summation over possible abstract triangulations. The UV cut-off is the 
length a of the sides of the simplices. Using this formalism, one can formulate a 
Euclidean theory of quantum gravity using as building blocks Euclidean equilateral 
simplices and one obtains a lattice version of two-dimensional quantum gravity. It 
can be solved analytically for finite a and agrees with a continuum quantization of 
two-dimensional Euclidean gravity (quantum Liouville theory) in the limit a — >• 0. 

However, in spacetime dimension larger than two this Euclidean lattice approach 
does not seem to have the desired continuum limit. This apparent failure was a key 
motivation for introducing a modified approach based on so-called causal dynamical 
triangulations (CDT). It realizes, in a nonperturbative context, ideas put forward 
in earlier work [4], which advocated that in a gravitational path integral with the 
correct, Lorentzian signature of spacetime one should sum over causal geometries 
only. 



2.2 The combinatorial solution in two dimensions 

Let us describe the explicit solution of CDT in two dimensions, that is, one space 
and one time dimension [7]. 

The model is defined as follows. The topology of the underlying manifold is 
taken to be S l x [0, 1], with "space" represented by the closed manifold S l . We con- 
sider the evolution of this space in "time". No topology change of space is allowed. 

The geometry of each spatial slice is uniquely characterized by the length as- 
signed to it. In the discretized version, the length L will be quantized in units of a 
lattice spacing a, i.e. L = I ■ a where / is an integer. A slice will thus be defined by 
I vertices and / links connecting them. To obtain a 2d geometry, we will evolve this 
spatial loop in discrete steps. This leads to a preferred notion of (discrete) "time" 
t , where each loop represents a slice of constant t . The propagation from time-slice 
t to time-slice t + 1 is governed by the following rule: each vertex i at time t is 
connected to kj vertices at time t + 1, kj > 1, by links which are assigned a "time- 
like" squared edge length —a 2 . The kj vertices, kj > 1, at time-slice f + 1 will be 
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connected by k, — 1 consecutive space-like links, thus forming k, — 1 triangles. Fi- 
nally the right boundary vertex in the set of kj vertices will be identified with the left 
boundary vertex of the set of kj + \ vertices. In this way we get a total of £J =1 (fc, ,— 1) 
vertices (and also links) at time-slice t + 1 and the two spatial slices are connected 
by L'=i h — h + k+i triangles, see Fig. 1. 




Fig. 1 The propagation of a spatial slice from time step / to step f + 1 in two-dimensional causal 
triangulations. The ends of the strip should be joined to form a band with topology S 1 x [0, 1]. 



The elementary building blocks of a geometry are therefore triangles with one 
space- and two time-like edges. We define them to be flat in the interior. A consistent 
way of assigning interior angles to such Minkowskian triangles is described in [5]. 
The angle between two time-like edges is y„ = — arccos|, and between a space- 
and a time-like edge y s , = § + i arccos §> summing up to y tt + 2y st = %. The sum 
over all angles around a vertex with j incoming and k outgoing time-like edges (by 
definition j,k > 1) is given by 2tt+ (4 — j — k) arccos |. The regular triangulation 
of flat Minkowski space corresponds to jf = k = 2 at all vertices. The volume of a 
single triangle is given by ^r-a . 

One may view these geometries as a subclass of all possible triangulations that 
allow for the introduction of a causal structure. Namely, if we think of all time-like 
links as being future-directed, a vertex v' lies in the future of a vertex v iff there is 
an oriented sequence of time-like links leading from v to v'. Two arbitrary vertices 
may or may not be causally related in this way. 

In quantum gravity one sums over all geometries connecting, say, two spatial 
boundaries of length L\ and Lq,, with the weight of each geometry g given by 



m 



S[g]=Ao U=g (in 2d), 



(1) 



where Aq is the bare cosmological constant. If in the discretized model we have a 
piecewise linear geometry made from n triangles the corresponding action will be 



S=A 



V5a 2 



Xn, 



:A, 



V5a 2 



(2) 



In our discretized model the boundaries will be characterized by integers l\ and 1%, 
the number of vertices or links at the two boundaries. The path integral amplitude 
for the "propagation" from geometry l\ to h will be the sum over all interpolating 
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surfaces of the kind described above, with a weight given by the discretized version 
of (1). Let us call the corresponding amplitude G^ (hjfa)- Thus we have 

G ( l\h,h) = iG { l ) {l u h;t), (3) 

t=\ 

G { l\h,h;t) = i,G { l\h,l;\)lG { l\l,l 2 ,t-\), (4) 

i=i 

G { l\hM;l)= l T E d^i* (5) 

where X denotes the bare dimensionless lattice cosmological constant 3 defined in 
(2), and where t denotes the total number of time-slices connecting l\ and 1%, 

From a combinatorial point of view it is convenient to mark a vertex on the 
entrance loop in order to get rid of the factors / and 1 // in (4) and (5), that is, 

G x (h,l 2 ;t)=hG^\h,h;t). (6) 

(The unmarking of a point may be thought of as the factoring out by (discrete) 
spatial diffeomorphisms). Note that G x (h >h\ 1) plays the role of a transfer matrix, 
satisfying 

G x {l u hA+t2)=Y. G dh,l\tx)G x {lMt2) (7) 

i 

G x (l h l 2 ;t+l)=^G x (k,l;l)G x (l,l 2 U)- (8) 

Knowing Gx(h : h'A) allows us to find Gx{h:h',t) by iterating (8) t times. This 
program is conveniently carried out by introducing the generating function for the 
numbers Gx(h,h'j), 

G x (x,y;t)=Y,x k y'Gx(k,l;t), (9) 

k.l 

which we can use to rewrite (7) as 

G l (x,y-t i +t 2 )= j ^-G x {x,z- X ;h)G x {z^ti), (10) 

where the contour should be chosen to include the singularities in the complex z- 
plane of G x {x, z~ ' ; 1 1 ) but not those of G x (z, y ; ?2 ) ■ 

One can either view the introduction of G x (x,y;t) as a purely technical device 
or take x and y as related to boundary cosmological constants. Let A, and Xf denote 
dimensionless lattice boundary cosmological constants, such that if the entrance 



3 One obtains the renormalized (continuum) cosmological constant A by an additive renormaliza- 
tion, see below. 



CDT — an Entropic Theory of Quantum Gravity 7 

boundary consists of k links the lattice boundary action will be X{k, or, introducing 
a dimensionful bare lattice boundary cosmological constant A, = A, i/a and a contin- 
uum boundary length L, = ka, A{Lj (and similarly Aq = Xq/ci etc.)- We now write 



A- = e' A <=e' A ' a , 



y = e' A /=e ,Aoa , (11) 



such that x k = e ik becomes the exponential of the boundary cosmological term, 
and similarly for y 1 — e f l . Let us for notational convenience define 

g = e x . (12) 

For the technical purpose of counting we view x,y and g as variables in the complex 
plane. In general the function 

G(x,y;g;t)=G x {x,y;t) (13) 

will be analytic in a neighbourhood of (x,y,g) = (0,0,0). 

From the definitions (5) and (6) it follows by standard techniques of generating 
functions that we may associate a factor g with each triangle, a factor x with each 
vertex on the entrance loop and a factor y with each vertex on the exit loop, leading 
to 

G (, y ; 8 ;i) = |(4fo)')'-|w'= (1 . gj) ^. gy) . (14) 

Formula (14) is simply a book-keeping device for all possible ways of evolving 
from an entrance loop of any length in one step to an exit loop of any length. The 
subtraction of the term 1/(1 — gx) has been performed to exclude the degenerate 
cases where either the entrance or the exit loop is of length zero. 
From (14) and eq. (10), with t\ = 1, we obtain 

G(x,y;g;t) = -^G(-?— ,y;g;t-l). (15) 

l-gx l-gx 

This equation can be iterated and the solution written as 

2 

G(x,y;g;t) = F?(x)F}(x)---F t 2 _ l (x)- * xy (16) 

[1 - gF, - 1 (*)] [1 - gF t - 1 (x) - gy] 

where F t (x) is defined iteratively by 

F 'W = i — l-rv F °W= X - ( 17 > 

l-gF t -i(x) 

Let F denote the fixed point of this iterative equation. By standard techniques one 
readily obtains 
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l-xF + F 2, -Hx-F) 1- Jl-Ag 2 

F ^= F i- xF+ F^l-Fr F= 2 8 - (18) 

Inserting (18) in eq. (16), we can write 

F 2 '(l~F 2 ) 2 xy 

G{x,y;g,t)= 7 - — — —. , — r (19) 

(A t - B,x) (A, - B, (x +y)+ C,xy) 

= F 2 '(l~F 2 ) 2 xy ^ (2Q) 

(l-xF)-F 2t + 1 (F-x) (l-xF)(l-yF)-F 2t (F-x)(F-y) 

where the time-dependent coefficients are given by 

A, = 1-F 2l+2 , B t =F{\-F 2, ) 1 C t =F 2 (l-F 2 '- 2 ). (21) 

The combined region of convergence to the expansion in powers g k x l y m , valid for 
all t is 

\g\<\, M<i, bl<i- (22) 

The asymmetry between x and y in the expressions (19) and (20) is due to the 
marking of the entrance loop. If we also mark the exit loop we have to multiply 
G x (h,h;t) by l 2 . We define 

Gf{l u h;t)=hG x {h,h;t) = hl 2 G { l\l u h;t). (23) 



The corresponding generating function G^ 2 '(x,y;g;t) is obtained from G(x,y;g;t) 
by acting withy ^, 

<*<***)- H^-J?" r (24) 

(A, -B,(x+y)+C,xyy 

We can compute G x {hih',t) from G(x,y;g;t) by a (discrete) inverse Laplace 
transformation 

G ^ l ^-f^xf^y^y^ G ^^ (25) 

where the contours should be chosen in the region where G(x,y;g;t) is analytic. A 
more straightforward method is to rewrite the right-hand side of (19) as a power 
series in x and y, yielding 

F^l-F 2 ) 2 ^ ™"(^2H l l+ l 2 -k-l f A,Q\ k 

G x[h,h;t) - h Al&l2+2 ^ Jt!(/,-Jt-l)!(/ 2 -jfc-l)! \^f) ' 

(26) 
which, as expected, is symmetric with respect to l\ and l 2 after division by l\. 
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In the next section we will give explicit expressions for Gx(h,h\t), Gx{h,h) 
and Gx (x,y) (the integral of Gx (x,y\ t) over t) in a certain continuum limit. 



2.3 The continuum limit 

The path integral formalism we are using here is very similar to the one used to re- 
present the free particle as a sum over paths. Also there one performs a summation 
over geometric objects (the paths), and the path integral itself serves as the propa- 
gator. From the particle case it is known that the bare mass undergoes an additive 
renormalization (even for the free particle), and that the bare propagator is subject 
to a wave-function renormalization (see [6] for a review). The same is true in two- 
dimensional gravity, treated in the formalism of dynamical triangulations [6]. The 
coupling constants with positive mass dimension, i.e. the cosmological constant and 
the boundary cosmological constants, undergo an additive renormalization, while 
the partition function itself (i.e. the Hartle-Hawking-like wave functions) undergoes 
a multiplicative wave-function renormalization. We therefore expect the bare cou- 
pling constants X , Xj and Ao to behave as 



4^- 2 = %+a, a,.^ = ^+1, a„^ = ^ 

4 or cr a a a a 



where A,X,Y denote the renormalized cosmological and boundary cosmological 
constants and where we have absorbed a factor v5/4 in the definition of A. 
If we introduce the notation 

g c = e lC *, Xc = e iC \ y c = e C V, (28) 

for critical values of the coupling constants, it follows from (11) and (12) that 

g = g c e? A , x = x c e iajt , y = y c e af . (29) 

The wave-function renormalization will appear as a multiplicative cut-off dependent 
factor in front of the bare "Green's function" G(x,y;g;t), 

G A (X,Y;T) = ]ima r >G(x,y;g;t), (30) 

a— >o 

where T = at, and where the critical exponent 77 should be chosen such that the 
right-hand side of eq. (30) exists. In general this will only be possible for particular 
choices of g c , x c and y c in (30). 

The basic relation (7) can survive the limit (30) only if 77 = 1, since we have 
assumed that the boundary lengths L\ and L2 have canonical dimensions and satisfy 
Lj = alj. 

From eqs. (19) and (21) it is clear that we can only obtain a nontrivial continuum 
limit if \F\ —} 1. This leads to a one-parameter family of possible choices 
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1 



2 cos a 



for F = e m , aeIR, (31) 



for critical values of g. It follows from (12) that most values of g c correspond to 
a complex bare cosmological constant A. However, the renormalized cosmological 
constant A in (27) (depending on how we approach g c in the complex plane) could 
in principle still be real. 

A closer analysis reveals that only at g c = ±1/2, corresponding to a = 0,71, is 
there any possibility of obtaining an interesting continuum limit. Note that these 
two values are the only ones which can be reached from a region of convergence of 
G(x,y;g;t) (see Fig. 2). Note also that requiring the bare A to lie inside the region 




Fig. 2 The circle of convergence in the complex g plane (radius 1/2), and the critical lines, ending 
in g = ±l/2. 



of convergence when g — > g c leads to a restriction ImA > on the renormalized 
cosmological constant A, since \g\ < A =>• Im A > ln2. 

Without loss of generality, we will consider the critical value g c = 1/2. It corre- 
sponds to a purely imaginary bare cosmological constant X c := Cxi a 1 = —i\n2/a 2 . 
If we want to approach this point from the region in the complex g-plane where 
G(x,y;g;t) converges it is natural to choose the renormalized coupling A imaginary 
as well, A = iA, i.e. 

X = i — f + /A. (32) 

a L 

One obtains a well-defined scaling limit (corresponding to A e IR) by letting A — ¥ 
k c along the imaginary axis. The Lorentzian form for the continuum propagator is 
obtained by an analytic continuation A — > —/A in the renormalized coupling of the 
resulting Euclidean expressions. 

At this stage it may seem that we are surreptitiously reverting to a fully Euclidean 
model. We could of course equivalently have conducted the entire discussion up to 
this point in the "Euclidean sector", by omitting the factor of —i in the exponential 
(1) of the action, choosing A positive real and taking all edge lengths equal to 1. 
However, from a purely Euclidean point of view there would not have been any 
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reason for restricting the state sum to a subclass of geometries admitting a causal 
structure. The associated preferred notion of a discrete time allows us to define an 
"analytic continuation in time" (we will discuss this in more detail later for higher- 
dimensional gravity). Because of the simple form of the action in two dimensions, 
the rotation 



dx dt\]-g\ or -> i / dx dteu^fgeu (33) 

to Euclidean metrics in our model is equivalent to the analytic continuation of the 
cosmological constant A . What is special about the above situation is that we per- 
form the analytic continuation configuration by configuration, i.e. geometry by ge- 
ometry. That is possible because of the particular set of causal geometries we have 
chosen to include in the regularized path integral. Moreover, as will be clear later, it 
is a feature which extends to higher dimensions too: each piecewise linear geometry 
with Lorentzian signature we use in the path integral has an analytic continuation 
to a Euclidean piecewise linear geometry and one has a relation like (33) for the 
Einstein-Hilbert actions of the two geometries (see later for details). 

From (19) or (20) it follows that we can only get macroscopic loops in the 
limit a — > if we simultaneously take x,y — > 1. (For g c = —1/2, one needs to 
take x,y —> — 1. The continuum expressions one obtains are identical to those for 
g c = 1/2.) Again the critical points correspond to purely imaginary bare bound- 
ary cosmological coupling constants. We will allow for such imaginary couplings 
and thus approach the critical point A, = A/ = from the region of convergence of 
G(x,y;g;t), i.e. via real, positive X,Y where 

A,- = iXa, A/ = iYa. (34) 

Again X and Y have an obvious interpretation as positive boundary cosmological 
constants in a Euclidean theory, which may be analytically continued to imaginary 
values to reach the Lorentzian sector. 
Summarizing, we have 



\z- Aal -> \{\ - \Aa\ (i.e. F = 1-aVX) 



(35) 



as well as 



e- Xa -> 1 - aX , y = <T aY -> 1 - aY, (36) 



where the arrows — > in (35) and (36) should be viewed as analytic coupling constant 
redefinitions of A,X and Y, which we have performed to get rid of factors of 1/2 
etc. in the formulas below. With the definitions (35) and (36) it is straightforward 
to perform the continuum limit of G(x,y;g,t) as (x,y,g) — > (x c ,y c ,g c ) = (1,1,1/2), 
yielding 



G A (X,Y;T) = 



\A e- 2 ^ 



(VX +x) +c- 2 ^ t (Va -x) 
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1 



(VA +X)(VA +Y) -e- 2 ^ T (VA -X)(VA -Y) 
For T — >• °° one finds 



(37) 



G A (X,Y;T) T -^ ^J- ^. (38) 

From Ga (X, Y; T) we can finally calculate G A {L\,L 2 ; T), the continuum ampli- 
tude for propagation from a loop of length L\ , with one marked point, at time-slice 
T = to a loop of length L 2 at time-slice T, by an inverse Laplace transformation, 

G A (L l ,L 2 \T)= f dX J dYe XLl e YL2 G A {X,Y;T). (39) 

J — ioQ J — ZOO 

This transformation can be viewed as the limit of (25) for a —¥ 0. The continuum 
version of (10) thus reads 

G A (X,Y;T 1 +T 2 )= f dZG A (X,-Z;Ti)G A (ZJ;T 2 ), (40) 

J — /co 

where it is understood that the complex contour of integration should be chosen to 
the left of singularities of G A (X , — Z; T\ ) , but to the right of those of G A (Z, Y, T 2 ) . 
By an inverse Laplace transformation we get in the limit T — >• °° 

G A {L U L 2 ;T) ^4L ie -^ L ' +L 2) e- 2 ^ T , (41) 

where the origin of the factor L\ is the marking of a point in the entrance loop. For 
r->0we obtain 

G A (X,Y;T) T -^ j±j, (42) 

in agreement with the expectation that the inverse Laplace transform should behave 
like 

G A {L U L2;T) T -^8(L X -L 2 ). (43) 

The general expression for G A (L\,L 2 ;T) can be computed as the inverse Laplace 
transform of formula (37), yielding 

e -[c°thv^T]v/A(L 1+ L 2 ) ^XIJ^ /2VAL^.\ 

sinh VAT L 2 VsinhVAT / ' 

where I\ (x) is a modified Bessel function of the first kind. The asymmetry between 
L\ and L 2 arises because the entrance loop has a marked point, whereas the exit 
loop has not. The amplitude with both loops marked is obtained by multiplying 
with L 2 , while the amplitude with no marked loops is obtained after dividing (44) by 
L\ . The highly nontrivial expression (44) agrees with the loop propagator obtained 
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from a bona-fide continuum calculation in proper-time gauge of pure 2d gravity by 
Nakayama [33]. 

The important point we want to emphasize here is that the additive renormaliza- 
tion of the cosmological constant is an entropic effect when calculated after rotation 
to Euclidean signature. In fact, we can write the propagator (13) as 

G(x,y,g;t)=^ X k y l g" E TTFv < 45 > 

k,l,n T{k,l,n) C( ^ > 

where the summation is over all causal triangulations T(k,l,n) (as defined above 
and rotated to Euclidean signature), consisting of n triangles and with the two 
boundaries made of k and I links. C(T) is the order of the so-called automorphism 
group of graph T and in our case, with a mark on one boundary, C(T) = 1. The 
critical point is g c = 1/2. That can only be the case because the number of (causal) 
triangulations constructed from n triangles grows exponentially as e" ln2 . The con- 
tinuum renormalized cosmological constant, as defined by eq. (35), emerges when 
taking the difference between the value of the action for a geometry made of n trian- 
gles and the entropy of the configurations with a given action (which in this case is 
proportional to the number of triangles n. More precisely, let the number of causal 
triangulations which can be constructed from n triangles be 

Jf(n) = /(«) e^", Xc = ln2, (46) 

where /(«) is a prefactor growing slower than exponentially, and which can also 
depend on the boundary cosmological constants x,y, a dependence we will suppress 
here. We can now write eq. (45) as 

G(A)=£/(»)e-( A - A > g = *-*■. (47) 

Introducing the notation A = na 2 for the continuum area (again disposing of a factor 
V5/4 for notational simplicity) we see that (35) can be written as 

A=A f +Afl 2 , (48) 

introducing the renormalized cosmological constant A . Eq. (47) can now be written 

as roo 

2\ „-AA 



G(A)= / 6Af(A/a') e - M , (49) 

Jo 

with the continuum action AA and the nontrivial physics contained in the function 
/(A/a 2 ). 

The two-dimensional CDT model can be generalized in a number of ways: one 
can use different weights and explore the universality of the model [8] and there 
exists a Hamiltonian formulation [9]. Matter can be coupled to the model [10] and 
one can weaken the causal constraints [11]. In addition, one can relax the constraint 
of a bounded geometry [12]. 
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3 Causal dynamical triangulations in four dimensions 
3.1 The choice of triangulations 

The generalization from two spacetime dimensions to three or four is in principle 
straightforward [13, 14]. In what follows we will concentrate on the 4d case. We 
consider spacetimes with the topology [0, 1] x S 3 . In principle we can choose any 
spatial topology, as long as we do not allow it to change during time evolution. Here, 
for simplicity, we will always take the topology of space to be that of a three-sphere. 

Suppose now that we have a foliation of spacetime where "time" is taken to 
mean proper time. Each time-slice, with the topology of S 3 , is represented by a 
three-dimensional triangulation. We choose as the set of possible triangulations of 
S 3 those which can be constructed from gluing together tetrahedra whose links are 
all of length a s = a, playing the role of lattice spacing and UV cut-off. These tetra- 
hedra are thus building blocks for our curved ^-geometries, which we take to be 
piecewise linear. The curvature of such a piecewise linear geometry is located at the 
links. A number of tetrahedra will share a link. Each tetrahedon has a dihedral angle 
associated with that link, and the sum of dihedral angles of the tetrahedra sharing 
the link would add up to 2n if space was flat around that link. If the dihedral angles 
add up to something different it signals that the piecewise linear space is not flat. 
In our case the tetrahedra are all identical with a dihedral angle f / = arccos(l/3), 
which implies there is no exact tessellation of flat three-dimensional space using 
equilateral tetrahedra. However, it is not important for the use we are making of the 
piecewise linear geometries: we use them in the path integral where we sum over 
all geometries (of a given, fixed topology). Thus the important question is whether 
the set of piecewise linear geometries we are using is dense in the set of geome- 
tries relevant for the path integral. To answer this question, we need to know the 
measure on the set of geometries. Presently we do not even have a mathematical 
characterization of the set of geometries to be used in the path integral (and the path 
integral has of course not been defined in any mathematical sense). Will the set of 
relevant geometries be like the set of paths used in the path integral of the particle, 
i.e. all continuous paths? Will the path integral over geometries include all "continu- 
ous" geometries? Because the answer is presently unknown, we will proceed with a 
straightforward generalization of the class of piecewise straight paths of the particle 
case and see what we get. 

We now connect two neighbouring ^-triangulations 73(1) and Ti(2), associated 
with two consecutive discrete proper times 1 and 2, and create a four-dimensional, 
piecewise linear geometry, such that the corresponding four-dimensional "slab" con- 
sists of four-simplices, has the topology of [0, 1] x 5 3 , and has 73(1) and 73(2) as 
its three-dimensional boundaries. The spatial links (and subsimplices) contained in 
these four-dimensional simplices lie in either 73(1) or 73(2), and the remaining links 
are time-like with proper length squared aj = —ad 1 , a > 0. Subsimplices which 
contain at least one time-like link we will call "time-like". In discrete units, we can 
say that 73(1) and 73(2) are separated by a single step in time direction, correspond- 
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(4,1) 



(3,2) 



Fig. 3 A (4, 1)- and a (3,2)-simplex interpolating between two neighbouring spatial slices. The 
time-reversed (1,4)- and (2, 3)-simplices are obtained by turning these upside down. 



ing to a time-like distance ^jda in the sense that each link in the slab which con- 
nects the two boundaries has a squared proper length —aa 2 . It does not imply that 
all points on the piecewise linear manifold defined by T$(l) have a proper distance 
squared —aa 2 to the piecewise linear manifold defined by T}(2) in the piecewise 
Minkowskian metric of the triangulation. 

Thus, our slabs or "sandwiches" are assembled from four-dimensional simplicial 
building blocks of four kinds, which are labelled according to the number of ver- 
tices they share with the two adjacent spatial slices of constant integer proper time. 
A (4, l)-simplex has one tetrahedron (and consequently four vertices) in common 
with 73(1), and only one vertex in common with T$(2). It has four time-like links, 
connecting each of the four vertices in 73 (1 ) to the vertex belonging to Tj(2). In the 
second kind of four-simplex, of type (1,4), the roles of 73(1) and 73(2) are inter- 
changed. By contrast, a (3,2)-simplex has a spatial triangle (and consequently three 
vertices) in common with the slice 73(1) and a spatial link (with two vertices) in 
common with 73 (2), together with six time-like links connecting the two slices. The 
corresponding (2,3)-simplex is again obtained by interchanging 73(1) and 73(2). 
The allowed simplices, up to time reversal, are shown in Fig. 3. For ourpurposes, we 
need not keep track of the numbers of four-simplices and their time-reversed coun- 
terparts separately, and will denote the total number of (4, 1)- and (l,4)-simplices 
by AQ and similarly the total number of (3,2)- and (2, 3 )-simplices by AQ ' . An 
allowed four-dimensional triangulation of the slab has topology [0, 1] x S 3 , is a sim- 
plicial manifold with boundary, and is constructed according to the recipe above. To 
summarize, a path in the gravitational path integral consists of a sequence of trian- 
gulations of S 3 , denoted by T$(k), k = 0,...,n, where the space between each pair 
73 (k) and Tj(k+l) has been filled in by a layer of four-simplices. In the path integral 
we sum over all possible sequences {73 (A:)} and all possible ways of triangulating 
the slabs in between Ti(k) and Tj(k+l). 
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3.2 The choice of action 

Piecewise linear geometries allow for a natural realization of the Einstein-Hilbert 
action, as discovered by Regge [15]. For a piecewise linear geometry in d dimen- 
sions, defined by a triangulation with length assignments to each link, the curvature 
is concentrated on the (d — 2)-dimensional sub-simplices. For example, in three di- 
mensions the curvature is located at the links, as described above. In dimension four, 
the curvature is concentrated at the triangles of the triangulation. A direct measure 
of the sectional curvature of the subspace perpendicular to a given triangle is the 
deficit angle one can associate with it, defined as the difference between the sum 
of "dihedral" angles of the four-simplices sharing the triangle and 2%. A deficit an- 
gle different from zero signifies nonvanishing curvature. In our case we have two 
kinds of triangles, purely space-like ones and time-like ones (where two of the links 
are time-like). The local contribution to the total integrated curvature from such a 
(d — 2) -dimensional sub-simplex is the volume of the sub-simplex multiplied by the 
deficit angle around it, and the integrated discrete (scalar) curvature action is then 
the sum of these contributions. This leads to a discretized Einstein-Hilbert action of 
the form 



l_|dV=^)(*(x)-2A) -► 
S Resse = k( £ Vo1(A)t(2w- £ ©)+ E Vol(A) (2tt- ^ 



S EH 

16nG. 



space— like 4— simplices time— like 4— simplices 

A at A A at A 

-A( £ Vol(4,l)+ £ Vol(3,2)), (50) 

(4,I)&(1.4)- (3,2)&(2.3)- 

tetrahedra tetrahedra 

where R is the Ricci scalar curvature and A the cosmological constant. Furthermore, 
we can read off that the constant k is proportional to the inverse of the gravitational 
coupling constant G and the constant A is proportional to A jG. 

It is straightforward to calculate the volumes of the spatial and time-like triangles 
as well as the four-volumes of the (4,1)- and (3,2)-simplices. This enables us to 

(4 1) 

express the discretized Einstein-Hilbert action as function of the numbers N\ ' and 
/y(3,2) defined in the previous subsection and the total number A^o of vertices in the 
triangulation, leading to 

5 Re gg e = (#„_£). knVAa + l + 



N 



f a) • ( k^\/4a+l-V3k arcsinh 



2 2\/2V3a+l 



3k n 2a + 1 . \/8a" 

- v 4a + 1 arccos — — A, - 



2 2(3a + l) 96 
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AQ ■ f kn\/4a+ 1 + ^-^- arcsinh 



2(3a + l) 



3£ -/ -1 4a + 3 

- V 4a + 12 arccos — - + arccos - 



4 V 2y/2y/2a + W3a + l 4(2a + l) 



For details of the derivation of this formula from (50) we refer to [ 14] . The quantity 
X denotes the Euler characteristic of the four-dimensional spacetime and appears be- 
cause we have been using the Euler relation No — N\ +N2 — N3+N4 = % along with 
other so-called Dehn-Somerville relations to express the numbers N\, Ni and N3 in 
terms of A/4 and Nq, where Ni counts the number of /-dimensional (sub-)simplices 
of a given triangulation. The constant a comes from allowing for a finite scaling 
aj = —aa 2 s between the length assignment of the proper length of space- and time- 
like links. 

The Lorentzian action (51) is obviously real for a > 0. We now want to study 
its rotation to Euclidean signature. This is naturally implemented by performing 
an analytic continuation in a from positive to negative a. In this way, the squared 
proper lengths of all links become positive and we have a piecewise linear geometry 
of Euclidean signature, where the links connecting two vertices from neighbouring 
time slices have length aj = \oc\a 2 . Of course, we also want the Euclideanized four- 
simplices to be nondegenerate. This requires — a > 7/12, the value below which 
a (3,2)-simplex becomes degenerate. Performing the analytic continuation in the 
complex lower a-half-plane and ending at a negative value smaller than -7/12 results 
in the following Euclidean action (note that we have made a redefinition a — >• —a, 
such that now (X > 7/12): 



sf^ e = ~knV^^T(N -x) 



,(4,1 



+W[ky/to=ll-~ 



% 



v/3 

: arcsin - 



2 V4a-1 2VV3a - 1 



3 2a-l\ „ V8a-3 
-arccos- +A- 



2 6a-2J 96 



\/3 6a -5 3 4a -3 

+Na kv4-a — 1 — % -A , arccos + - arccos 

4V4a - 1 6a-2 4 8a-4 



3 1 N , ^ ^MoTl \ 

-arccos — - + A . (52) 



2 2 v / 2V2a- V3a-1/ 96 

This action is precisely the Regge action for a piecewise linear geometry, con- 
structed from Euclidean building blocks where aj = \oc\a$. By analytic continuation 
in the complex lower-half a-plane we have therefore arrived at the usual formula 
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^Minkowski ^ ~^Euclidean (P3) 

from quantum field theory, where Seudidean expressed in terms of (continuum) Eu- 
clidean geometry is 

In quantum field theory on a flat, Minkowskian background such a "Wick rotation" 
does not usually refer to a specific configuration in the path integral. A configu- 
ration will typically not even be differentiable and therefore not have any analytic 
continuation. Rather, it will refer to the functional form of the action S[(j>] and how 
it changes by formally replacing itu by f£. This reflects the fact that the Feynman 
propagator, i.e. the result of a functional integration over the field(s) <p, has specific 
analytic properties which allow for a rotation in t to the corresponding Euclidean 
propagator. In our framework the situation is different, in that each individual piece- 
wise linear geometry can be rotated to a corresponding Euclidean geometry and its 
associated action transforms in accordance with (53). 

For what follows it will be convenient to simplify this action and write it as 

Sf sge = -(K + 6A)No + K 4 (N i 4 4A) +M 3 ' 2) ) +A(2N { 4 4A) +A/f 2) ), (55) 

where Kb = kK\J\u — 1, K4 is a linear combination of k and A with coefficients 
depending of a, and A is, for fixed Kb and K4, a function of a, chosen such that 
for a = 1 (i.e. when a, = a s ) we have 4=0. For the values of Kb and K4 which 
will be of interest for us any A > corresponds to a, < a s . In Fig. 4 we show a as 
function of A for Kb = 2.2 and K4 chosen critical (a concept to be discussed shortly). 
Furthermore, we have dropped the term involving %. It plays no role in the dynamics 
since we are not changing the topology. 



3.3 The entropic theory of gravity and its phase diagram 

Having chosen the cut-off and an action, we can now write down the path integral 
or partition function for the CDT version of quantum gravity, 

Z(G,A)=|f%]e- 5 I H k] -> Z{ko^A)=Y,^^ Se{T \ (56) 

where the summation is over all causal triangulations T of the kind described above, 
and we have dropped the superscript "Regge" on the discretized action Se given by 
(55). Like in the two-dimensional model, the factor 1/Ct is a symmetry factor, 
given by the inverse of the order Ct of the automorphism group of the triangulation 
T. Note that we can write the partition function as 
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Fig. 4 Plot of the asymmetry factor a, defined as aj = aa 2 s , where a, and a s are the lengths of 
time-like and space-like links, plotted as a function of A. The horizontal line is a = 7/12, the 
lowest allowed value of a, where the (3,2)-simplices collapse in the time direction. 



Z(K ,K 4 ,A)= £ e -(^)iV4 S -AN?» S (K,+6A)N £ 



n 4 A 4X 



.A/ 



r(«,^ 4,1) ,ivb) 



1 



(57) 



Introducing 



we can write 



x = e- (K *+ A \ y = e- A , z = e (K ° +6A \ 



(58) 



Z(x,y,z) = U *»< /4 ' ^o jk(N a ,N^' 1} ,N ) 



N 4 ,N^ [) ,N Q 



(59) 



where ,yV(N4,N\ ' ,Nq) denotes the number of CDT configurations with A/4 four- 

simplices of which AQ ' are of type (4, 1 ) or ( 1 , 4), and with /Vo vertices, including 
symmetry factors. Thus the calculation of the partition function is in principle a 
combinatorial problem, just as in two dimensions where we could solve the prob- 
lem explicitly (and formula (59) is similar to (45)). This is the reason why we call 
the model entropic: the partition function is entirely determined by the number of 
geometries in the simplest possible way, namely, by being the generating function 
for these numbers. The counting of geometric "microscopic" configurations of given 
(N4,N\ ' ,Nq) is their entropy in the sense of statistical models. Unlike in two di- 
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mensions, it has until now not been possible to solve this counting problem analyti- 
cally, which means that we will have to rely on numerical methods. 

Let us first understand better the nature of the partition function given by (57). 
We can write the sum as 

Z(k ,k 4i A) = £ e -(K*+^ Z N4 (Ko,A), (60) 

where Z# 4 (jq),4) is the partition function for a fixed number N4 of four-simplices, 
namely, 

z N4 (K Q ,A) = Y -L e -^V*) e e*+w(^ 4 ). (6i) 

t n , Ct 

One can show that Z^ 4 (Kb, A) is exponentially bounded as a function of N 4 [16], 

Z N4 (k q ,A) < e ^ +A ^f(N 4 ,K ,A), (62) 

where f(N 4 ) grows slower than exponentially. We call K^ the critical value of K4. 
It is a function of A and Ko and plays the same role as X c in the two-dimensional 
model discussed above: the partition function is only defined for K4 > fcj and the 
"infinite-volume" limit, where (N 4 ) — > °° can only be achieved for K\ — » k:|. We 
are interested in sending the lattice spacings a = a Sl a t to zero while keeping the 
physical four-volume, which is roughly N 4 a 4 , fixed. Thus we want to consider the 
limit N 4 — >■ °°, and fine-tune K4 to kJ for fixed Ko, A. This fine-tuning is similar to 
the fine-tuning X — >• X c in the two-dimensional model. Like there, we expect the 
physical cosmological constant A to be defined by the approach to the critical point 
according to 

* 4 ^ + 16^ 4 < (63) 

an equation similar to (48). It ensures that the term 

(K4-K c 4 )N 4 = -^—V 4 , V 4 =N 4 a\ (64) 

167TG 

gives rise to the standard cosmological term in the Einstein-Hilbert action. The cor- 
responding phase diagram is described qualitatively by Fig. 5. The shaded surface 
is the "critical surface", which we want to approach from above (by decreasing K4). 
We put "critical surface" in quotation marks since more correctly fine-tuning to this 
surface corresponds to taking the limit of infinite four- volume, which does not nec- 
essarily imply also a continuum limit. The situation here may be different from that 
of the two-dimensional model, where approaching X c automatically meant taking 
a continuum limit. Two-dimensional quantum gravity is of course a very simple 
model with no propagating degrees of freedom, whereas in four-dimensional quan- 
tum gravity we expect to have genuinely propagating field degrees of freedom. Thus 
the situation is more like in ordinary Euclidean lattice field theory /critical phenom- 
ena, where "infinite volume" does not necessarily mean "continuum limit". 
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Kf (K n ,A) 




Fig. 5 The phases A, B and C in the coupling constant space spanned by (kq,A, K4). Phase C is 
the one where extended four-dimensional geometries emerge. 



A good example of what one might expect is the Ising model on a finite lattice. 
To obtain a phase transition for this model one has to take the lattice volume to 
infinity, since there are no genuine phase transitions for finite systems. However, 
just taking the lattice volume to infinity is not sufficient to ensure critical behaviour 
of the Ising model. We also have to tune the coupling constant to its critical value, 
at which point the spin-spin correlation length diverges. Similarly, having placed 
ourselves on the "critical" surface of CDT quantum gravity or, rather, its "infinite- 
volume" surface, we can discuss the various phases, indicated in the figure as A, B 
and C. We move between these phases by changing the bare coupling constants A 
and K"o. The solid lines drawn on the surface are the phase transition lines between 
the (potentially) different phases, which can be of first or higher order. In order to 
go from the discrete lattice to the continuum theory, we are often interested in a 
second- or higher-order phase transition, since such transitions are associated with 
divergent correlation lengths of the field propagators, allowing one to forget about 
the lattice spacing relative to the correlation length. We will be looking for such 
transition points. 

How can one imagine obtaining an interesting continuum behaviour as a function 
of KqI For the purpose of illustration, let us assume that the subleading correction 
/(Af4,fco) has the form 
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f(N 4 ,K ) = e k ^^ r \ (65) 

(where later we will check numerically that such a term is indeed present, cf. eq. 
(78)). The partition function now has the form 

Z(K4, Kb) = ^e-^-^+^oJv^ (66) 

For dimensional reasons we expect the classical Einstein term in the action to scale 
like 



U*SVg(Z)m « #> (67) 



\bnGJ * V6V ^ vw G ' 

motivating the search for a value Kg with k(K^) = 0, with the approach to this point 
governed by 

*(*b)«^> i.e. k(K )VN^c^-. (68) 

With such a choice we can identify a continuum limit where (A/4), calculated from 
(66) (by a trivial saddle point calculation), goes to infinity while a —>■ 0, 

V*4) „ _l'<r._v*>W.J-H'v-. - l./Wr ~ /!/„ „<n? ** A?^4 - (""^ 



3-(*4-^)JV 4 +*(xb)v^J 4(K4-K£) 2 A 2 fl 4 ' 



Thus we find 



^4' z( ^^ )wex K4(?^))- exp (^)' (70) 

as one would naively expect from Einstein's equations, with the partition function 
being dominated by a typical instanton contribution, for a suitable constant c. 

The actual set-up for the computer simulations is slightly different from the the- 
oretical framework discussed above, in that we choose to work with a fixed num- 
ber of four-simplices N4 in the computer simulations, rather than fine-tuning K4 to 
its critical value. We can perform computer simulations for various N4 (and fixed 
Kq,A) and in this way check scaling with respect to A/4. This is an alternative to 
fine-tuning K4, and much more convenient from a computational point of view. For 
large A/4 we can then check whether there are any finite-size effects or whether ef- 
fectively we already are at the "critical surface" shown in Fig. 5. In addition, we 
fix the total number N t of spatial slices, with proper-time labels t\, t% = t\ + a t , up 
to % =h + (N t — l)a t , where At = a t is the discrete lattice spacing in the temporal 
direction 4 , and denote by T = N t a t the total extension of the universe in proper time. 
For convenience we identify fjy+i with t\, in this way imposing the topology S l x S 3 
rather than [0, 1] x S 3 . This choice does not affect physical results, as will become 
clear below when we present the numerical results. 



4 The separation between adjacent slices is a,, in the sense that all links connecting two neighbour- 
ing slices have length a,, as discussed above. 
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Fig. 6 The phase diagram of four-dimensional quantum gravity, defined in terms of causal dynam- 
ical triangulations, parametrized by the inverse bare gravitational coupling Kq and the asymmetry 
parameter A . 



Finally, the computer simulations are so-called Monte Carlo simulations, where 
the computer simply generates a sequence of configurations, in this case piecewise 
linear geometries, with the correct probability distribution, as dictated by the mea- 
sure and action used in the path integral. It is based on a local updating in terms of 
"moves", which change the piecewise linear geometry in a well-defined way (see 
[14] for a detailed description). These geometries are used to calculate expectation 
values of observables. 



3.4 The actual phase diagram 



Based on computer simulations with N4 = 80.000 we have constructed the phase 
diagram shown in Fig. 6 [23]. The dotted lines in the figure represent mere extra- 
polations, and lie in a region of coupling constant space which is difficult to access 
due to the inefficiencies of our computer algorithms. 

There are three phases, labeled A, B and C. In phase C, which had our main 
interest in [18, 19, 21, 20], we observe a genuinely four-dimensional universe in 
the sense that as a function of the continuum four- volume V4 (linearly related to the 
number of four-simplices), the time extent scales as V 4 and the spatial volume as 
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3/4 

V 4 . Moving into phase A, these scaling relations break down. Instead, we observe 
a number of small universes arranged along the time direction like "pearls on a 
string", if somewhat uneven in size. Individual beads along the time direction can 
grow and shrink, be created or disappear as a function of the Monte Carlo time 
used in the simulations. These small universes are connected by thin "necks", i.e. 
slices of constant integer time t n , where the spatial 5 3 -universes are at or close to the 
smallest three-volume permitted (consisting of five tetrahedra glued together), thus 
preventing "time" from becoming disconnected. 

By contrast, phase B is characterized by the "vanishing" of the time direction, in 
the sense that only one spatial hypersurface has a three-volume appreciably larger 
than the minimal cut-off size of five just mentioned. One might be tempted to con- 
clude that the resulting universe is three-dimensional, just lacking the time direction 
of the extended universe found in phase C. However, the situation is more involved; 
although we have a large three-volume collected at a single spatial hypersurface, 
the corresponding spatial universe has almost no extension in the spatial directions. 
This follows from the fact (ascertained through direct measurement) that it is pos- 
sible to get in just a few steps from one tetrahedron to any other by moving along 
the centres of neighbouring tetrahedra or, alternatively, from one vertex to any other 
along a chain of links. The Hausdorff dimension is therefore quite high, and possi- 
bly infinite. Let us assume for the moment that it is indeed infinite; then the universe 
in phase B has neither time nor spatial extension, and there is no geometry in any 
classical sense. 

We can now give the following qualitative characterization of the three phases 
in terms of what we will provisionally call "average geometry". The universe of 
phase C exhibits a classical four-dimensional background geometry on large scales, 
such that (geometry) 7^ 0. One may even argue that (geometry) = const, in view 
of the fact that according to the minisuperspace analysis of [20, 21, 36] and al- 
lowing for a finite, global rescaling of the renormalized proper time, the universe 
can be identified with the round S 4 , a maximally symmetric de Sitter space of con- 
stant scalar curvature (as we will describe in detail below). By contrast, in phase 
B the universe presumably has no extension or trace of classicality, corresponding 
to (geometry) = 0. Lastly, in phase A, the geometry of the universe appears to be 
"oscillating" in the time direction. The different behaviour of typical configurations 
is shown in Fig. 7. The "time" direction is horizontal and we plot the three-volume 
Ni(t), i.e. the number of tetrahedra in a given time slice, as the circumference of a 
circle. The three phases are separated by three phase transition lines which meet in 
a triple point as illustrated in Fig. 6. 
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Fig. 7 The volume profiles of typical configurations in the phases A, B and C. Phase C (bottom 
figure) is the one where extended four-dimensional geometries emerge. 

4 The macroscopic de Sitter universe (phase C) 
4.1 Identifying the infrared part of the universe 



Phase C in the above-mentioned phase diagram has our main interest, because this 
is where we observe an extended four-dimensional universe. We will now discuss in 
more detail the geometric properties of this "macroscopic" universe. The measure- 
ments reported in this Section have been performed at the values {kq,A) = (2.2, 0.6) 
of the bare coupling constants, a point which lies well inside phase C. 

The Monte Carlo simulations referred to above will generate a sequence of space- 
time histories. An individual spacetime history is not an observable, in the same way 
as a path x(t) of a particle in the quantum-mechanical path integral is not. However, 
it is perfectly legitimate to talk about the expectation value (x(t)) as well as the 
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fluctuations around (x(t)), both of which are in principle calculable in quantum 
mechanics. 

Obviously, there are many more dynamical variables in quantum gravity than 
there are in the particle case. We can still imitate the quantum-mechanical situation 
by picking out a particular one, for example, the spatial three-volume V3 (t) at proper 
time t. We can measure both its expectation value (V3(f)) as well as fluctuations 
around it. The former gives us information about the large-scale "shape" of the 
universe we have created in the computer. A "measurement" of Vi(t) consists of a 
table A3 (/), where i = l,...,N t and A, denotes the total number of time slices. The 
time axis has a total length of N t time steps, where N t = 80 in the actual simulations, 
and we have cyclically identified time slice N t + 1 with time slice 1. 

What we observe in the simulations is that for the range of discrete volumes A4 
under study the universe does not extend (i.e. has appreciable three-volume) over 
the entire time axis, but rather is localized in a region much shorter than 80 time 
steps. Outside this region the spatial extension Nj(i) will be minimal, consisting 
of the minimal number (five) of tetrahedra needed to form a three-sphere S 3 , plus 
occasionally a few more tetrahedra. 5 This thin "stalk" therefore carries little four- 
volume, which means that in a given simulation we can for most practical purposes 
consider the total four- volume of the remainder, the extended universe, as fixed. 

In order to perform a meaningful average over geometries which explicitly refers 
to the extended part of the universe, we have to remove the translational zero mode 
present, see [21] for a discussion of the procedure. Having defined the "centre of 
volume" along the time direction of our spacetime configurations, we can now per- 
form superpositions of configurations and define the average (Ni(i)) as a function 
of the discrete time i. The results of measuring this average discrete spatial size of 
the universe at various discrete times i are illustrated in Fig. 8 and can be succinctly 
summarized by the formula 

^ BWll, ^i^™ , b)> s °* a59 ' <71) 

where A3 (/) denotes the number of three-simplices in the spatial slice at discretized 
time i and A4 the total number of four-simplices in the entire universe, % is a constant 
referring to the fact that we have a nonvanishing asymmetry A, which implies differ- 
ent lengths for time- and space-like links and consequently different four-volumes 
for four-simplices of type (4,1) and (3,2). Likewise, the ratio AQ ' /AQ ' depends 
on the choice of bare coupling constants and has to be measured. Of course, formula 
(7 1 ) is only valid in the extended part of the universe where the spatial three- volumes 
are larger than the minimal cut-off size. 

The data shown in Fig. 8 have been collected at the couplings (kq,A) = (2.2,0.6) 
and for A4 = 362.000. For these particular values of (kq,A) we have verified rela- 
tion (71) for A 4 ranging from 45.500 to 362.000 building blocks (45.500, 91.000, 



5 This kinematic constraint ensures that the triangulation remains a simplicial manifold in which, 
for example, two rf-simplices are not allowed to have more than one (d — l)-simplex in common. 
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Fig. 8 Background geometry (Ni ((')): MC measurements for fixed N4 = 362.000 and best fit 
(71) yield indistinguishable curves at given plot resolution. The bars indicate the average size of 
quantum fluctuations. 



181.000 and 362.000). After rescaling the time and volume variables by suitable 
powers of N4 according to relation (71), and plotting them in the same way as in 
Fig. 8, one finds almost total agreement between the curves for different spacetime 
volumes, as illustrated in Fig. 9. This constitutes a beautiful example of finite-size 
scaling. At least with regard to measuring the average three-volume Vi(t) all our 
discretized volumes N4 are sufficiently large to be treated as infinite, in the sense 
that no further changes will occur for larger N4. 

By contrast, the quantum fluctuations indicated in Fig. 8 as vertical bars for each 
discrete time i are volume-dependent and will become (relatively) larger when the 
the total four- volume is decreased. Eq. (71) shows that spatial volumes scale ac- 
cording to /V 4 ' and time intervals according to Ay , as one would expect for a 
genuinely /oMr-dimensional spacetime. This is exactly the scaling we have used in 
Fig. 9. It strongly suggests a translation of (71) to a continuum notation. The most 
natural identification is given by 



v^n'( t) ^v 4 ^Q 



(72) 



where we have made the identifications 



B 



SqK 



1/4 : 



Ati^V 3 (ti) = 2C 4 N 3 (i)a 4 , 



(73) 



such that 



dty/teV 3 (t) = V 4 . 



(74) 
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Fig. 9 Rescaling of time and volume variables according to relation (71) for A/4 = 45.500, 91.000, 
181.000 and 

N3(i)/NT- 



181.000 and 362.000. The plot also include the curve (71). More precisely: cr <* i/N 1 ^ and P(a) 



In (73), yfgTt is the constant proportionality factor between the time t and genuine 
continuum proper time T, T = ^fgT t t. (The combination Atiy/gnVj, contains Q, re- 
lated to the four- volume of a four-simplex rather than the three-volume correspond- 
ing to a tetrahedron, because its time integral must equal V4). Writing V 4 = %K 2 R A /2>, 
and y/gH = 7?/B, eq. (72) is seen to describe a Euclidean c/e Sitter universe (a four- 
sphere, the maximally symmetric space for positive cosmological constant) as our 
searched-for, dynamically generated background geometry! In the parametrization 
of (72) this is the classical solution to the action 



5 = 



1 



24nG 



dt 



V 3 (t) 



-k 2 vl'\t) 



■AV 3 (0 



(75) 



where k 2 = 9(27T 2 ) 2 ' 3 and A is a Lagrange multiplier, fixed by requiring that the 
total four-volume be V4, Jdt^/g^ Vj(t) = V4. Up to an overall sign, this is precisely 
the Einstein-Hilbert action for the scale factor a (f ) of a homogeneous, isotropic uni- 
verse (rewritten in terms of the spatial three-volume Vj(t) — 2n 2 a(t) 3 ), although 
we of course never put any such simplifying symmetry assumptions into the CDT 
model. The intriguing possibility of describing the data in terms of the minisuper- 
space model (75) was first reported in [36]. 

A discretized, dimensionless version of (75) is 
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^f (N 3 (i+l)-N 3 (i)) 2 r 1/3 N 

■W = fc l L ( jy-?^ + ^3 (0 ] > ( 76 > 

where £2 x &2- This can be seen by applying the scaling (71), namely, N?,(i) = 
N 4 ni(sj) and Sj = i/N,' . This enables us to finally conclude that the identifi- 
cations (73) when used in the action (76) lead naively to the continuum expression 
(75) under the identification 6 

a 2 \/Ca 5k 
G = 7-^r4r- W 

We also note that the reported scaling (71) implies that we can write 






\ m{Si+l) -^ {Si) )\hny\ Si ) 



1/4 

where ds, = l/N 4 . Thus, referring to (65), we see that there are at least terms in 
log /(Mt, kq,A) which scale like y/N^. To obtain the form (65) we need additional 
terms of entropic nature since the coefficient of y/TU, in (65) is assumed positive. 



4.2 The size of the universe and the flow of G 

It is natural to identify the coupling constant G multiplying the effective action 
for the scale factor with the gravitational coupling constant G. The effective action 
describing our computer-generated data is given by eq. (75), and its dimensionless 
lattice version by (76). The computer data allow us to extract k\ « a 2 /G, with a the 
spatial lattice spacing, and the precise constant of proportionality given by eq. (77). 
For the bare coupling constants (kq,A) = (2.2,0.6) we have high-statistics mea- 
surements for N4 ranging between 45.500 and 362.000 four-simplices (equivalently, 
AQ ' ranging between 20.000 and 160.000 four-simplices). The choice of A deter- 
mines the asymmetry parameter a, and the choice of (kq,A) determines the ratio £, 

(3 2) (4 1) 

between N 4 ~ ' and N\ . This in turn determines the "effective" four- volume C\ of 
an average four-simplex, which also appears in (77). The number So in (77) is de- 
termined directly from the time extension r un ; v of the extended universe according 
to 

? (4,l)\'/4 



Tuniy = 7CSo[N^ J j . (79) 

Finally, from our measurements we have determined k\ = 0.038. Taking everything 
together according to (77), we obtain G rj 0.23a 2 , or ipi sa 0.48a, where ipi = \/G 
is the Planck length. 



6 Due to the difference in four- volume between N A ' and N\ for a 7^ 1 we have to introduce a 
compensating factor so = so (N4 ) ' ' 4 / (AQ ' ) ' ' 4 , see [21] for a detailed discussion. 
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From the identification of the volume of the four-sphere, V4 = &7t 2 R 4 /3 = 
C4AQ ' a 4 , we obtain that R = 3. la. In other words, the linear size KR of the quan- 
tum de Sitter universes studied here lies in the range of 12-21 Planck lengths for 
N4 in the range mentioned above and for the bare coupling constants chosen as 
(Ko,A) = (2.2,0.6). 

Our dynamically generated universes are therefore not very big, and the quan- 
tum fluctuations around their average shape are large, as is apparent from Fig. 8. 
The presence of such fluctuations is evident in the bottom snapshot picture of the 
extended universe shown in Fig. 7, whose volume profile deviates from that of a 
regular sphere. The point is of course that we have to perform an averaging process 
to obtain the expectation value of the volume profile, and this is precisely what we 
have been doing numerically. It is rather surprising that the semiclassical minisuper- 
space formulation gives an adequate description - at least for the volume profile - 
for universes of such a small size, a fact that should be welcome news to anyone per- 
forming semiclassical calculations to describe the behaviour of the early universe. 
However, when looking at more local geometric properties of the universe, our lat- 
tices are still coarse compared to the Planck scale lp\ because the latter corresponds 
to roughly half a lattice spacing. If we are after a theory of quantum gravity valid on 
all scales, we are specifically interested in uncovering phenomena associated with 
Planck-scale physics. In order to collect data which are free from unphysical short- 
distance lattice artefacts at this scale, we would ideally like to work with a lattice 
spacing much smaller than the Planck length, while still being able to set by hand 
the physical volume of the universe studied on the computer. The way to achieve 
this is by changing the bare coupling constants Kq,A such that the coefficient K in 
(77), G =: Ka~, is changed to a larger value. However, A" is a combination of a num- 
ber of factors and they might change differently when Kb, A are changed. It is thus 
a (computer-)experimental exercise to find a path in the (kq,A) coupling-constant 
plane such that K increases and we end up with Ipi = \f~Ka ^> a. We will discuss 
later whether such a path exists. 



5 Constructive evidence for the effective action 

We have found a perfect fit (7 1) to the emergent background geometry and the curve 
can be related to the continuum effective action (75). However, it is still of interest 
to investigate to what extent the action (76) can be rederived from the data. Interest- 
ingly, as we shall see below, this can largely be done. 

The data at our disposal are: (i) the measurement of the three-volume A3 (/) at the 
discrete time step i, and of the three-volume correlator NifyN^tf). Having created 
Q statistically independent configurations A3 (/) by Monte Carlo simulation allows 
us to construct the average 
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iV 3 (0:=(iV3(0> = 7^1^(0, (80) 

9=1 



where the superscript in (•p q > denotes the result of the q'th configuration sampled; 
(ii) the covariance matrix 

C(i,j) = ^Z,(4 q) (i)-N 3 (i)Wi q) U)-N3U)y (81) 

We now assume we have a discretized action which can be expanded around the 
expectation value Nj, (?) according to 

S di scr[N + n] = S discr [N] + \Y. n ' P 'i n J + <9 (" 3 )' < 82 > 

If the quadratic approximation describes the quantum fluctuations around the expec- 
tation value N well, the inverse of the operator P will be a good approximation to 
the covariance matrix. Conversely, still assuming the quadratic approximation gives 
a good description of the fluctuations, the P constructed from the covariance matrix 
will to a good approximation allow us to reconstruct the action via (82). 

Fig. 10 shows the measured covariance matrix C(i, j) and its inverse, the operator 
P. Some care is needed in inverting C(i,j) since it has two zero modes, one from 
the constraint that N4 is kept fixed, and (an approximate) one from the fact that 
the translational mode of the 'centre of volume' can only be fixed up to a lattice 
spacing, see [21] for a detailed discussion. As is clear from the figure, the inverse 
P is completely dominated by the stalk data. This feature is unavoidable: while the 
correlation matrix is dominated by long-range fluctuations, the inverse matrix will 
be dominated by short-distance fluctuations, i.e. the fluctuations in the stalk, which 
by definition are associated with cut-off energies. 

Looking at the inverse P of the measured covariance matrix, we observe that 
to very good approximation it is small and constant, except on the diagonal and 
the entries neighbouring the diagonal. This means that we can decompose it into a 
"kinetic" and a "potential" term. The kinetic part P km is defined as the matrix with 
non-zero elements on the diagonal and in the neighbouring entries, such that the 
sum of the elements in a row or column is always zero. The potential part P pot is 
then given by whatever remains along the diagonal. We therefore arrive at a tentative 
representation of P as 

P iJ =fff n +P%°<, (83) 

Pfj n =PiAij, Plf'^mdij, (84) 

where the matrices Ay and 8{j are essentially defined through the construction just 
described 7 . We know P from the data, and can make a least-^ 2 fit to determine the 



7 For details of normalization and subtleties in the definition of P k "' and P'"" related to the zero 
modes we refer to [21]. 
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Fig. 10 The covariance matrix C (top) and its inverse (bottom). 



numbers p, and «,. For details we refer again to [21]. The results are shown in Figs. 
11 and 12. 

Let us look at the discretized minisuperspace action (76) which has served as in- 
spiration for the definition of P kln and P p< " . Expanding N^(i) to second order around 
Nj,(i), one obtains the identifications 



N 3 (i) = —> U"(N 3 (i)) 
Pi 



-Uj, 



(85) 
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Fig. 11 The directly measured expectation values N^(i) compared to the averages N}(i) recon- 
structed from (85), for Kq = 2.2 and A = 0.6. 



1/3 

where U(Ni(i)) = kikaNJ (i) denotes the potential term in (76). We use the fitted 
coefficients pi to reconstruct^ (/) and then compare these reconstructed values with 
the averages Nj,(i) measured directly. Similarly, we can use the measured m,'s to re- 

- — S/3 

construct the second derivatives U"(Ni (/)) and compare them to the form N 3 (/) 
coming from (76). 

The reconstruction of Nj(i) is illustrated in Fig. 11 for a given four- volume N4 
and compared with the directly measured expectation values N^(i). One observes 
that the reconstruction works very well and, most importantly, that the coupling con- 
stant k\, which in this way can be determined independently for each four- volume 
N4, really is independent of N4 in the range of Afy's we have considered, as it should 
be. 

We will now try to extract the potential U" (JV3 (/) ) from the information contained 
in the matrix P pot . The determination of U"(N^(i)) is not an easy task as can be 
understood from Fig. 12, which shows the measured coefficients m,- extracted from 
the matrix P pot , and which we consider rather remarkable. The interpolated curve 
makes an abrupt jump by two orders of magnitude going from the extended part of 
the universe (stretching over roughly 40 time steps) to the stalk. The occurrence of 
this jump is entirely dynamical, since no distinction has ever been made by hand 
between stalk and bulk. In order to extract physical information related to a genuine 
potential like the one appearing in (76), we of course must restrict ourselves to 
the region inside the "blob", corresponding to the data range encircled in Fig. 12. 
From the figure it is also clear that extracting U"(Ni(i)) from the data available is a 
nontrivial task. 
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160.000. Data chosen from the encircled region are independent of short- 



The range of the discrete three-volumes A3 (i) in the extended universe is from 
several thousand down to five, the kinematically allowed minimum. However, the 
behaviour for the very small values of A3 (/) near the edge of the extended universe 
is likely to be mixed in with discretization effects. In order to test whether one 
really has a Ay (/)-term in the action one should therefore only use values of N^(i) 
somewhat larger than five (shown as the encircled region in Fig. 12). This has been 
done in Fig. 13, where we have converted the coefficients m, from functions of the 
discrete time steps i into functions of the background spatial three-volume A3 (/) 
using the identification in (85) (the conversion factor can be read off the relevant 
curve in Fig. 11). The data presented in Fig. 13 were taken at a discrete volume 
A4 = 160.000, and fit well the form A 3 , corresponding to a potential ^2^3 . 

-5/3 

Apart from obtaining the correct power A 3 for the potential for a given space- 
time volume A4, it is equally important that the coefficient in front of this term be 
independent of A/4. This seems to be the case as is shown in Fig. 14, where we have 
plotted the measured potentials in terms of reduced, dimensionless variables which 
make the comparison between measurements for different A4's easier. - In sum- 
mary, we conclude that the data allow us to reconstruct the action (76) with good 
precision. 
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6 Connection to Hof ava-Lifshitz gravity 



While we have verified that the action (76) describes the data well, one aspect of 
the formula may lend itself to a more general interpretation. Until now we have 
chosen to view the "experimental" formula (71) as describing a round four-sphere 
via the identifications (72) and (73). However, the potential asymmetry between 
space and time introduced in our model by working with a time foliation allows for 
a different interpretation, namely, that space and time really behave differently, as 
we will explain in what follows. 

Although at the level of the piecewise linear structures we have a precise rela- 
tion between the coupling constant A and the asymmetry parameter a, this relation 
enters the construction only in a relatively weak way, in that for a given A the bare 
action we use is the Regge-Einstein-Hilbert action for a piecewise linear manifold 
with the given connectivity and the length assignment af = aa 2 s . Nowhere else does 
this length assignment appear. Of course, ;/ the model had a well-defined pertur- 
bative expansion, one could have chosen the bare coupling constants such that the 
classical action (and, by implication, the relation between a s and a t ) played an im- 
portant role in the path integral. However, as already explained in the introduction, 
this is not the case. Rather, our choice of bare coupling constants is dictated by the 



36 



J. Ambj0rn, A. Gorlich, J. Jurkiewicz and R. Loll 



18 
16 
14 
12 
10 



l 



Linear Fit 

Monte Carlo, 160k 

Monte Carlo, 80k 

Monte Carlo, 40k 




Fig. 14 The dimensionless second derivative u 
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plotted against v 5 / 3 , where v = 

= 40.000, 80.000 and 160.000, Kb = 2.2 
and A = 0.6. One expects a universal straight line near the origin (i.e. for large volumes) if the 
power law U(Ns) °= N 1 ^ is correct. 



N3/N4 is the dimensionless spatial volume, for AQ ' 



wish to find nontrivial physics, which restricts us to a region far away from where 
the action term dominates the entropy of configurations. 

In other words, the effective action we have obtained bears only indirect traces of 
the classical action put into the path integral. Similarly, the precise relation between 
space and time directions in the final continuum theory will be determined by the 
statistical averages resulting from the full path integral, rather than the parameter a 
put into the discretized action. 

To illustrate the nonperturbative mechanism at work, let us consider the measure- 
ments in phase C, where we observe a macroscopic universe whose extension in the 

1/4 

time direction scales like iV 4 ' for different, fixed four-volumes N4. It should be 
emphasized that this scaling behaviour is by no means predetermined, for example, 
time in phase B scales completely differently (in fact, the time extension vanishes 
there, which is of course an extreme situation). One is therefore led to conclude that 
time and (the linear extension of) space scale identically in phase C. This is corrob- 
orated by evidence that well inside phase C there seems to be a well-defined notion 
of a "physical" proper time extent independent of a, as one would have expected 
naively [22]. 
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However, the fact that the behaviour in phase B is very different leaves open 
the interesting possibility that a nontrivial scaling relation between time and spa- 
tial extent may ensue when one approaches the B-C phase transition line. In that 
case one could try to describe the situation by a continuum effective action where 
time and space have different dimensionality. This is precisely what P. Hofava at- 
tempted in his novel class of gravity theories [23], a topic we will return to shortly. 
But even in phase C, where time and space scale in the same way, one could ex- 
plore the consequences of relaxing the notion that time and space should be related 
exactly as they are in general relativity. What we have done until now is to inter- 
pret our results in the continuum limit in terms of the classical theory. This can be 
achieved by making a global rescaling of continuum proper time, which we have di- 
rect access to through the preferred time foliation. More specifically, we have been 
checking our data against the cosmological minisuperspace model (75), which can 
be derived from general relativity by assuming spatial homogeneity and isotropy. 
Eq. (75), when written in terms of the scale factor 8 a(x), T denoting proper time, 
reads 

dr (flfl 2 + a-Afl 3 ). (86) 



4 

All our measurements are perfectly consistent with an effective action of this form, 
but the "cosmological" observables we have been considering so far cannot dis- 
criminate between this and more general cosmologies coming from a generalized 
"gravity" theory with a built-in anisotropy between space and time, like Hofava's. 
In the latter, one can again assume spatial homogeneity and isotropy to obtain cos- 
mological solutions, which in the Euclidean sector, in the infrared limit, arise from 
an action quite similar to (86) [24, 25, 26], namely, 

5=| /dTf3(3A-l)flfl 2 + 67fl-fl 3 (6A+y(a))V (87) 

The "potential" V(a) in this expression has an expansion in inverse powers of a, 
coming from the higher-order spatial derivative terms in the Hofava-Lifshitz action. 
In the actual computer simulations, in both (86) and (87) A is a Lagrange multiplier 
rather than a cosmological constant, which ensures that the four-volume is kept 
fixed. As long as we can fix proper time only up to a constant and as long as we 
cannot measure reliably the correction term V(a) containing inverse powers of a, it 
is not really possible to distinguish "experimentally" between (86) and (87) in terms 
of a reconstruction of the action, as we did in the previous Section. In this situation, 
whenever A > 1/3 and y > 0, a rescaling of time in (87) leads to the same form 
as (86) up to a constant of proportionality. Taking into account the difficulties in 
verifying the mere existence of the linear term in (86) from the data (cf. Fig. 12), it 
is clear that we cannot presently extract from the data in a reliable way a potential 
V(a) that depends on inverse powers of a, starting with a~ 4 . The only region in 
Fig. 12 where such information could be extracted is for the smallest values of a. 



8 equivalently, one can also work with the associated three-volume V(t) = 2^ 2 a(r) 3 as the basic 
configuration space variable 
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Unfortunately, this is also where lattice artefacts will be important and artefacts from 
the "stalk" may get mixed in with genuine continuum physics; a glance at Fig. 12 
reveals how large the effects of the stalk are right next to the small-a region. 

This not withstanding, we can discuss the qualitative correspondence between 
the Hofava scenario and our phase diagram. In our earlier analysis of the different 
phases of CDT quantum gravity, we have chosen for a particular qualitative descrip- 
tion to match precisely that of a Lifshitz phase diagram [27, 28]. The qualitative 
feature we want to emphasize in this context is that the role played by "average ge- 
ometry" in quantum gravity bears an intriguing resemblance to that played by the 
Lifshitz field 0. In an effective Lifshitz theory, the Landau free energy density F(x) 
as function of an order parameter (x) takes the form 9 



F(x)=a 2 §(xy+a^{xy+a 6 4>{xf + ...+c 2 {d a $y+d 2 {d fi $) z +e 2 { 

(88) 
where for a rf-dimensional system a = m + 1 , . . . , d, j3 = 1 , . . . , m. Distinguishing 
between "a"- and "j3 "-directions allows one to take anisotropic behaviour into ac- 
count. For a usual system, m — and a phase transition can occur when a 2 passes 
through zero (say, as a function of temperature). For «2 > we have <j> = 0, while for 
«2 < we have 1 1 > (always assuming 04 > 0). However, one also has a transition 
when anisotropy is present (m > 0) and d 2 passes through zero. For negative d 2 one 
can then have an oscillating behaviour of in the m "j3"-directions. Depending on 
the sign of a 2 , the transition to this so-called modulated or helical phase can occur 
either from the phase where = 0, or from the phase where 1 1 > 0. We conclude 
that the phases C, B, and A of CDT quantum gravity depicted in Fig. 8 can be put 
into one-to-one correspondence with the ferromagnetic, paramagnetic and helical 
phases of the Lifshitz phase diagram 10 if we in place of we use "average geom- 
etry". The triple point where the three phases meet is the so-called Lifshitz point, 
where in the Lifshitz model one can have a nontrivial scaling. 

The critical dimension beyond which the mean-field Lifshitz theory alluded to 
above is believed to be valid is d c =4 + m/2. In lower dimensions, the fluctuations 
play an important role and so does the number of components of the field 0. This 
does not necessarily affect the general structure of the phase diagram, but can alter 
the order of the transitions. Without entering into the details of the rather complex 
general situation, let us just mention that for m = 1 fluctuations will often turn the 
transition along the A-C phase boundary into a first-order transition. Likewise, most 
often the transition between phases B and C is of second order. 

We conclude that the structure of the Lifshitz phase diagram is presently com- 
patible with our CDT observations. Based on the order parameter investigated in 
[17], the A-C transition of CDT quantum gravity looks like a clear-cut first-order 
transition. On the other hand, the verdict is still out for the B-C line. The signals 
from the Monte Carlo simulations are ambiguous, and it appears that close to the 



9 see, for example, [28] for an introduction to the content and scope of "Landau theory" 

10 For definiteness, we are using here a "magnetic" language for the Lifshitz diagram. However, 
the Lifshitz diagram can also describe a variety of other systems, for instance, liquid crystals. 
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transition line our algorithms for updating the geometries need to be improved to 
produce fully reliable results. 

Phase transitions of higher than first order are of intrinsic interest since they may 
serve as points where one can define the continuum theory. Since it appears that we 
have a well-defined infrared limit in phase C, such points would naturally be UV 
fixed points, and moving away from them should bring us to the IR limit. Given 
the overall structure of the phase diagram, the possible scenarios are as follows: if 
either the A-C or the B-C line is second-order, we can potentially use any point on it 
to attempt to define a UV limit. By contrast, if they are first-order lines, we are left 
with two interesting points which may be associated with a higher-order transition: 
the endpointfb of the open B-C line (cf. Fig. 6), where the phase transition may be 
of higher order than along the line itself, and the Lifshitz triple point P t , where the 
transition may also be of higher order. 

Assuming there were such higher-order phase transition points, how would we 
determine whether the U V limit of the theory is isotropic in space and time or, more 
generally, of Hofava-type? One defining aspect of Hofava-Lifshitz gravity is the 
assumption that the scaling dimensions of space and time differ in the ultraviolet 
regime. This difference is used to construct a theory containing higher-order spatial 
derivatives in such a way that it is renormalizable. How would one observe such a 
difference in the present lattice approach? - Consider a universe of time extent T, 
spatial extension L and total four- volume Va{T,L). By measuring T and L we can 
establish the mutual relations 

T « vl' d \ L « (V 4 '- 1M ) ' M oc r(*-i)M. (89) 

Well inside phase C we have measured d, = 4 and d s = 3, in agreement with what is 
expected for an ordinary four-dimensional spacetime. If the dimension [T] of time 
was z times the dimension [L] of length, we would have 

(90) 



dt-l 



We observed earlier that well inside phase B both d s and d t must be large, if not 
infinite. In case the B-C phase transition is second-order, it may happen that z goes 
to a value different from 1 when we approach the transition line. To investigate 
this possibility, we have tried to determine z as a function of the parameter A as 
A — >• 0. For A > 0.3 one obtains convincingly d t « 4 and d s ~ 3 and thus z ~ 1. We 
can make an even stronger statement, namely, that the data does not contradict the 
interpretation of A as an (unphysical) asymmetry parameter: when A is increased 
the corresponding a is decreasing (see Fig. 4), while the number of lattice spac- 
ings in the time direction is increasing, which at least qualitatively allows for the 
interpretation that the physical "time" is independent of a. By contrast, for A < 0.3 
the quality of our results does not allow for any definite statements. Autocorrelation 
times become very long and there may be large finite-volume effects, which obscure 
the measurements and which are precisely based on finite-size scaling. 
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To summarize, there is still a distinct possibility that a nontrivial scaling with 
z 7^ 1 can occur when one approaches the transition line between phases B and 
C. Should the B-C transition turn out to be a first-order transition, the interesting 
points would be Pq and P t . It would then be natural to conjecture that the point P t is 
a Hofava-Lifshitz point with a nontrivial scaling relation between space and time. 
The point Pq, which does not appear in a standard Lifshitz diagram, would perhaps 
be a natural candidate for an isotropic scaling point. On the other hand, if the B- 
C transition line is second-order, it leaves open the interesting possibility that the 
critical exponent changes continuously from z = 3 at the Lifshitz point P t to z = 1 at 
the (hypothetically) isotropic point Pq. 



7 Making contact with asymptotic safety 

As we discussed earlier, it is presently difficult to get close to the B-C phase tran- 
sition line, which is needed if we want to achieve a higher resolution in the UV 
regime, such that the lattice spacing is much smaller than the Planck length. Also, 
we do not know yet whether in such a limit we have isotropy in space and time, 
like in the asymptotic safety approach, or need to invoke an anisotropic scenario 
as outlined above. For the time being, let us assume that the endpoint Pq of the B- 
C transition line in the phase diagram of Fig. 6 corresponds to an isotropic phase 
transition point. How can one make contact with the gravitational renormalization 
group treatment? The standard way would be to "measure" observables (by lattice 
Monte Carlo simulations), like a mass in QCD or the string tension in Yang-Mills 
theory. For definiteness, consider the string tension, which has mass dimension two. 
The measurements, for some choice go of the bare coupling constant, will give us a 
number G (go). We now write 

o{go) = <?Ra 2 (go), (91) 

where Gr is the physical string tension and a (go) describes the dependence of the 
lattice spacing a on the bare coupling constant. Being able to write down a relation 
like this for all observables, where a (go) is determined by the renormalization group 
equation 

a^ = -P(g ), (92) 

da 

allows us to define a continuum theory at a fixed point g^ where B (go) = 0, since 
there we can take a (go) — > when go — > gp. In the case of QCD or Yang-Mills theory 
the fixed point is the Gaussian fixed point g^ = 0, but in the more general setting of 
asymptotic safety it will be non-Gaussian. 

Assume now that we have a fixed point for gravity. The gravitational coupling 
constant is dimensionful, and we can write for the bare coupling constant 

G(a)=a 2 G(a) 1 a— = -B(G), B(G) = 2G-cG 3 + • • • . (93) 

da 
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The IR fixed point G = corresponds to G constant while the putative non-Gaussian 
fixed point corresponds to G — » G*, i.e. G(a) -» G*a 2 . In our case it is tempting to 
identify our dimensionless constant k\ with 1 jG, up to the constant of proportion- 
ality given in (77). Close to the UV fixed point we have 

G(a) = G*-Ka E , ki=k\+Ka s , c = -j3'(G*). (94) 

Usually one relates the lattice spacing near the fixed point to the bare coupling con- 
stants with the help of some correlation length <jj . Assume that ^ diverges according 

4to>) = ra <95) 

in the limit as we tune the bare coupling constant go —>■ g^. This correlation length 
is associated with a field correlator and usually some physical mass m p /, by means 
of 

—jr- — ^- =m ph (a\ni -n 2 \) =m ph \x\ -x 2 \, (96) 

where \n\ — n 2 \ is a discrete lattice distance and \x\ — x 2 \ a physical distance. Re- 
quiring the physical quantities \x\ —x 2 \ and m„/, to remain constant as a — >• then 
fixes a as a function of the bare coupling constant, 

a = — \go-g* Q \ v - (97) 

cm ph 

Eq. (97) is only valid close to the fixed point and should be compared to the renor- 
malization group equation (92), from which we deduce that v = — l/|/3'(gg)|- 

In the gravitational case at hand we do not (yet) have observables which would 
allow us to define meaningful correlation lengths. At any rate, it is by no means 
a settled issue how to define such a concept in a theory where one integrates over 
all geometries, and where the length is itself a function of geometry (see [32] for 
related discussions). Instead, we construct from our computer-generated "data" an 
effective action, where all degrees of freedom, apart from the scale factor, have been 
integrated out. We impose the constraint that the data are fitted to a universe of total 
lattice four-volume A/4. Measurements are performed at different, fixed values of A/4, 
all the while maintaining the relation 11 

V 4 =N 4 a 4 . (98) 

We then "remove the regulator", by investigating the limit A/4 — » 00. In ordinary 
lattice field theory, we have two options for changing A/4; either we keep a fixed, 
and therefore change V4, or we keep V4 fixed and change a. Let us illustrate the 
difference in terms of a scalar field on a lattice. Its dimensionless action can be 
written as 



11 in principle, we should be taking into account the different volumes of the two types of four- 
simplices, which depend on A , but we will ignore these details to streamline the presentation 
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s=L(5>(*+M)-*(0) 2 +«<to 2 (o)> (99) 

i I* 

where i labels discrete lattice points and /x unit vectors in the different lattice direc- 
tions. The correlation length is approximately 1 /mo lattice spacings. Holding a fixed 
and increasing N4 is not going to change the correlation length in any significant way 
if N4 is sufficiently large. Thus the interpretation for fixed a is straightforward: the 
physical volume V4 is simply increased and finite-size effects will become smaller 
and smaller. However, we can also insist on an interpretation where V4 is kept fixed, 
N4 is increased and a decreased accordingly. In this case, the lattice becomes finer 
and finer with increasing N4. But now the physical interpretation of (99) will change 
with increasing N4, even if no bare coupling constant is changed, and the correlation 
length is still approximately 1/ma lattice spacings. Since the physical lattice length 

1/4 

a decreases proportional to 1/Ay the physical correlation length is going to zero, 
and the physical mass to infinity. This can be made explicit in (99) by introducing 
the lattice spacing a, 



(100) 



The physical mass is m„f, = mo/a and goes to infinity unless we adjust niQ. The 
factor 1 J a in front of the sum can be reabsorbed in a redefinition of if desired. 

In our case it is natural to consider V4 as fixed if we want to make contact with 
the continuum framework of asymptotic safety, since this will allow us to vary a. 
Suppose we have identified a fixed point which we consider interesting, e.g. the 
point Pq in our phase diagram. We now approach this point and measure k\. If it is 
a UV fixed point, eq. (94) tells us what to expect. Using (98), we can now convert 
this into an equation involving N4, and suitable for application in CDT simulations, 

h{N 4 )=k c l -KN; E/ \ (101) 

When we measured k\ (N4) deep inside phase C (at the point (K2,A) = (2.2,0.6)), 
we did not find any A^-dependence of k\ . However, according to the insights just 
presented, we should observe such a dependence at or close to a UV fixed point. 
As already noted earlier, an explicit verification of such a relation will have to await 
more reliable computer simulations close to the phase transition lines. 

In fact, we have already seen indications in CDT quantum gravity of a short- 
distance behaviour like that occurring in the asymptotic safety scenario. Recall the 
"naive" renormalization conditions (64) and (68). They were introduced mainly to 
illustrate how a renormalization procedure could lead to finite renormalized cos- 
mological and gravitational constants, both with a semiclassical interpretation. If 
we are close to the UV fixed point, we know that G will not be constant when we 
change scale, but G will. Writing G(a) = a 2 G*, eqs. (64) and (68) are changed to 

K4-K^^a 2 , *i(K§) = ^. (102) 
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Fig. 15 The data points along the central curve show the spectral dimension D$ ( o" ) of the universe 
as function of the diffusion time o". Superimposed is a best fit, the continuous curve Ds(o) = 
4.02— 1 19/(54+0"). The two outer curves quantify the error bars, which increase linearly with c 
(Measurements taken for a quantum universe with 181.000 four-simplices.) 



The first of these relations now looks two-dimensional (cf. eq. (48))! Nevertheless, 
the expectation value of the four- volume still satisfies the correct relation 



<V 4 ) = <^4)« 4 -^, 



(103) 



as follows from (69). 

Further hints of a two-dimensional signature at short distances have come from 
measuring the so-called spectral dimension. Essentially, this is the dimension a dif- 
fusing "liquid" would experience in a spacetime with a typical quantum geometry of 
the kind appearing in the gravitational path integral, see the original article [29] for 
details. Fig. 15 recalls the result of measuring the spectral dimension, which appears 
to change nontrivially as a function of diffusion time. Since short diffusion times 
probe short distances, we can read off from the fit indicated that the short-distance 
spectral dimension is close to two. What is most intriguing is that the short-distance 
result Ds = 2 has since been found in two other quantum field-theoretic approaches, 
namely, the renormalization group treatment a la Reuter [30] and Hofava-Lifshitz 
gravity [31]. At the same time, it means that this particular observable cannot dis- 
tinguish between isotropic and anisotropic quantum gravity theories. 
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8 Discussion 

The CDT model of quantum gravity, which we have described in these lectures is 
extraordinarily simple. It implements the path integral over causal geometries with 
a global time foliation. In order to perform the integration explicitly, we introduced 
a grid of piecewise linear geometries, similar to how one proceeds when defining 
the path integral in ordinary quantum mechanics. Next, we rotated each of these 
geometries to Euclidean signature and used as our bare action the Einstein-Hilbert 
action in Regge form. In terms of ingredients, that's all. 

The resulting superposition exhibits a nontrivial scaling behaviour as function of 
the four- volume, and we observe the appearance of a well-defined average geometry, 
that of de Sitter space, the maximally symmetric solution to the classical Einstein 
equations in the presence of a positive cosmological constant. The measurements 
performed so far definitely probe the quantum regime, since the fluctuations of the 
three-volume around de Sitter space are sizable, as can be seen in Fig. 8. Both the 
average volume profile and the quantum fluctuations around it are well described 
in terms of the minisuperspace action (75). A key feature to appreciate here is that, 
unlike in standard (quantum-)cosmological treatments, this description is the out- 
come of a nonperturbative evaluation of the full path integral, with everything but 
the scale factor (equivalently, Vo,(t)) summed over. Measuring the correlations of 
the quantum fluctuations in the computer simulations for a particular choice of bare 
coupling constants enabled us to determine the continuum gravitational coupling 
constant GasG« 0.42a 2 , thereby introducing an absolute physical length scale 
into the dimensionless lattice setting. Within measuring accuracy, our de Sitter uni- 
verses (with volumes lying in the range of 6.000-47.000 £%) are seen to behave 
perfectly semiclassically with regard to their large-scale properties. 

These semiclassical results "emerge" even if, as emphasized above, we are far 
from a region in coupling-constant space where the classical action can be consid- 
ered as dominant in the path integral. The resulting quantum geometry comes about 
through the interplay of the weight provided by the exponential of the bare action 
and the weight provided by the entropy of a particular kind of configuration. From 
this point of view the results are truly nonperturbative, even if they bear some simi- 
larity to the semiclassical minisuperspace results. It should also be emphasized that 
we have only derived an effective action for the scale factor, not for the "real", trans- 
verse degrees of freedom. The issue of how to do this for the latter remains to be 
addressed. 

The results we have reported are mostly "infrared" in nature. Our dynamically 
generated universes are macroscopic (although small) and - with the exception of 
the spectral dimension measurements - we are not yet probing Planckian (and possi- 
bly sub-Planckian) scales. It is a major issue whether such a short-distance comple- 
tion of gravity exists in a conventional, field-theoretical sense. "Asymptotic safety" 
is an attempt to formulate the general conditions for this to be the case. In its sim- 
plest realization on the lattice it requires a UV fixed point. This is precisely the kind 
of situation the CDT framework allows us to address; we have a phase diagram 
and potential fixed points, and as outlined above it is in principle possible to check 
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whether the UV behaviour is in accordance with the predictions from asymptotic 
safety. Critical slowing down close to the potential fixed points, i.e. a vast increase 
in the computer time needed to generate statistically independent geometries, has 
so far made it impossible to obtain reliable results, but work is now in progress to 
improve the updating algorithms. 

Once we succeed in approaching the phase transition line B-C in a reliable fash- 
ion, we may also be able to check whether the Hofava-Lifshitz scenario is realized 
in our model. In addition to the standard asymptotic safety picture, the CDT con- 
struction has the potential to accomodate also this scenario. Its obvious similarities 
with the defining properties of Hofava's "anisotropic gravity" approach are the dis- 
tinguished role of a time foliation, and the presence of a unitary time evolution (in 
CDT linked to the reflection positivity of the transfer matrix). 

In addition to this, we have already pointed to the striking similarity between 
the CDT and Lifshitz phase diagrams upon (somewhat loosely) identifying the Lif- 
shitz mean-field order parameter with "average geometry". If more specifically 
we want to relate the Lifshitz field to a mode of the gravitational field, the confor- 
mal factor appears as a natural candidate. The conformal mode is already known 
to play a decisive role in a variety of geometro-dynamical contexts. In the bare Eu- 
clidean Einstein-Hilbert action the kinetic term associated with the conformal factor 
appears with the wrong, negative sign, leading to ill-defined expressions for naive 
cosmological, Euclidean path integrals. In Euclidean noncritical string theory the 
dynamics of the conformal factor is believed to cause a transition from a "healthy 
phase" (where c < 1 for the central charge of matter) to a degenerate phase of so- 
called branched polymers (for c > l) 12 . A similar phenomenon was observed in the 
old Euclidean dynamical triangulations approach to quantum gravity, where the bare 
(inverse) gravitational coupling Kq plays the role of the central charge c, in the sense 
that for large values of Kq the configurations degenerate into branched polymers 
[39]. 

We have interpreted phase A (realized for large values of Kq) as the CDT remnant 
of the branched-polymer phase, likewise caused by the dominance of the conformal 
mode. This suggests that the A-C phase transition may be interpreted as a transition 
where the kinetic term of the conformal mode changes sign. This is precisely what 
happens at the A-C transition in a Lifshitz diagram, corresponding to the expression 
(88) for the free energy of the order parameter (p. 

The effective action for the conformal mode coming out of a nonperturbative 
gravitational path integral receives potential contributions from several sources: (i) 
from the bare action (where the kinetic conformal term has the "wrong", negative 
sign), (ii) from the measure, and (iii) from integrating out other field components 
and, where applicable, other matter fields. It has been argued previously that the 
Faddeev-Popov determinants obtained from gauge-fixing the gravitational path in- 
tegral contribute effectively with the opposite, positive sign to the conformal kinetic 
term [40, 41]. For example, when working in proper-time gauge, to imitate the time- 



12 In noncritical string theory there exists an analytic proof that for the dimension of (Euclidean) 
target spacetime d > 2 the string surfaces degenerate into branched polymers [34], see also [35] 
for similar results on a hypercubic lattice. 
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slicing of CDT, Euclidean metrics can be decomposed according to 13 

ds 2 = dr 2 + e 2 ^^gij(z, x )djdx j , (104) 

giving rise to a term — \ / G^e?^ y/Setg{d z <^) 2 in the bare gravity Lagrangian den- 
sity, where G^ b ' is the bare Newton's constant. According to [41], one expects that 
the leading contribution from the associated Faddeev-Popov determinant has the 
same functional form, but with a plus instead of a minus sign, and with a different 
dependence on G^ b ' . The presence of contributions of opposite sign to the effective 
action for the conformal mode (t,x) can therefore lead to two different behaviours, 
depending on the value of G^ b ' (equivalently, the Kq in our model), and thus account 
for the observed behaviour at the transition between phases A and C. 14 

In addition to the issues raised above, there is one more question which we would 
like to understand in more detail at this stage, which concerns the relation of our ef- 
fective gravitational coupling constant G to a more conventional gravitational cou- 
pling constant, defined directly in terms of coupling gravity to matter. It would be 
desirable to verify that defining the physical Newton's constant G as the coupling 
constant multiplying the effective action for the three-volume, as we have been do- 
ing so far, agrees with a gravitational constant defined more directly through matter 
coupling. In principle it is easy to couple matter to CDT quantum gravity, as we al- 
ready know from multiple studies in the Euclidean case [43], where spin, scalar and 
gauge fields have been considered. It is less straightforward to come up with a rea- 
sonably simple set-up for extracting the semiclassical effect of gravity on the matter 
sector (or vice versa), which is both well-defined on the ensemble of geometries 
and allows for effective computer measurements. Attempts in this direction were 
already undertaken in the "old" Euclidean approach [44, 45], and it is possible that 
similar ideas can also be used in our causal version of the theory. As a first step in 
this direction, the expected effect of a single point mass coupled to CDT quantum 
gravity on the volume profile of the universe has been quantified in [46]. Further 
work on coupled systems of matter and geometry is in progress. 
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13 The conformal decomposition of the spatial three-metric is essentially unique if one requires g/j 
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